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Single linear polymer chains in dilute solutions under good solvent conditions are studied by Monte 
Carlo simulations with the pruned-enriched Rosenbluth method up to the chain length N ~ 0(1O 4 ). 
Based on the standard simple cubic lattice model (SCLM) with fixed bond length and the bond 
fluctuation model (BFM) with bond lengths in a range between 2 and vTO, we investigate the 
conformations of polymer chains described by self-avoiding walks (SAWs) on the simple cubic lattice, 
and by random walks (RWs) and non-reversible random walks (NRRWs) in the absence of excluded 
volume (EV) interactions. In addition to flexible chains, we also extend our study to semiflexible 
chains for different stiffness controlled by a bending potential. The persistence lengths of chains 
extracted from the orientational correlations are estimated for all cases. We show that chains 
based on the BFM are more flexible than those based on the SCLM for a fixed bending energy. 
The microscopic differences between these two lattice models are discussed and the theoretical 
predictions of scaling laws given in the literature are checked and verified. Our simulations clarify 
that a different mapping ratio between the coarse-grained models and the atomistically realistic 
description of polymers is required in a coarse-graining approach due to the different crossovers to 
the asymptotic behavior. 


I. INTRODUCTION 

In the theoretical study of polymer physics 00 , com¬ 
puter simulations provide a powerful method to mimic 
the behavior of polymers covering the range from atomic 
to coarse-grained scales depending on the problems one 
is interested in 01. The generic scaling properties of 
single linear and branched polymers in the bulk or con¬ 
finement under various solvent conditions have been de¬ 
scribed quite well by simple coarse-grained lattice models 
(i.e., random walks (RWs), non-reversible random walks 
(NRRWs), self-avoiding random walks (SAWs), or inter¬ 
acting self-avoiding random walks (ISAWs) on a regular 
lattice, regarding the interactions between non-bonded 
monomers). As an alternative one can use coarse-grained 
models in the continuum, such as a bead-spring model 
(BSM) (where all beads interact with a truncated and 
shifted Lennard-Jones (LJ) potential while the bonded 
interactions are captured by a finitely extensible non¬ 
linear elastic (FENE) potential) using Monte Carlo and 
molecular dynamics simulations 0. On the one hand, 
however, as the size and complexity of a system increases, 
detailed information at the atomic scale may be lost when 
employing low resolution coarse-graining representations. 
On the other hand, the cost of computing time may 
be too high if the system is described at high resolu¬ 
tion. Therefore, more scientific effort has been devoted 
to developing an appropriate coarse-grained model which 
can reproduce the global thermodynamic properties and 
the local mechanical and chemical properties such as 
the intermolecular forces between polymer chains 00. 
While these models are already known since a long time, 
the present work is the first study presenting precise 
data on conformational properties of these model, when 
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a bond angle potential is included. 

In this paper we deal with linear polymer chains in 
dilute solutions under good solvent conditions, and de¬ 
scribe them by lattice models on the simple cubic lattice. 
Although coarse-grained lattice models neglect the chem¬ 
ical detail of a specific polymer chain and only keep chain 
connectivity (topology) and excluded volume, the uni¬ 
versal behavior of polymers still remains the same in the 
thermodynamic limit (as the chain length N -A oo) 0, 
Two coarse-grained lattice models, the standard simple 
cubic lattice model (SCLM) and the bond fluctuation 
model (BFM) [3. HI [71 . are considered for our simula¬ 
tions. The SCLM is often used for the test of new sim¬ 
ulation algorithms, and the verification of theoretically 
predicted scaling laws due to its simplicity and computa¬ 
tional efficiency. The BFM has the advantages that the 
computational efficiency of lattice models is kept and the 
behavior of polymers in a continuum space can be de¬ 
scribed approximately. The model thus introduces some 
local conformational flexibility while retaining the com¬ 
putational efficiency of lattice models for implementing 
excluded volume interactions by enforcing a single occu¬ 
pation of each lattice vertex. 

redThe excluded volume effect plays an essential role 
in any real polymer chain, while in a dilute solution un¬ 
der a theta solvent condition, or in concentrated polymer 
solutions such as melts, and glasses, the real polymer 
chain behaves like an ideal chain. The excluded volume 
constraint can easily be incorporated in the lattice mod¬ 
els by simply forbidding any two effective monomers oc¬ 
cupying the same lattice site (cell). Tries et al. have 
successfully mapped linear polyethylene in the melt onto 
the BFM and their results are in good agreement with 
experimental viscosimetric results quantitatively with¬ 
out adjusting any extra parameters El 0. Varying 
the backbone length and side chain length of the bottle¬ 
brush polymer based on the BFM, a direct comparison 
of the structure factors between the experimental data 
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for the synthetic bottle-brush polymer consisting of hy- 
droxyethyl methacrylate (PMMA) as the backbone poly¬ 
mer and flexible poly(n-butyl acrylate) (PnBA) as side 
chains in a good solvent (toluene) and the Monte Carlo 
results is given G2- Furthermore, the lattice models have 
also been used widely to inve stig ate the conformational 
properties of protein-folding [181 ] and DNA in chromo¬ 
somes fl9l - [2l[ in biopolymers. For alkane-like chains the 
angles between subsequent effective bonds are not contin¬ 
uously distributed, but only discrete angles are allowed. 
Therefore, the lattice model such as the SCLM where 
only the discrete angles 0° and 90° are allowed is an 
ideal model for studying alkane-like chains of different 
stiffnesses. 

A direct comparison of simulation results between the 
lattice models and the off-lattice BSM is also possi¬ 
ble, e.g. linear polymers [ 22 j and ring polymers (23| in 
a melt, adsorption of multi-block and random copoly¬ 
mers 24|, semiflexible chains under a good solvent condi¬ 
tion |25| , the crossover from semiflexible polymer brushes 
towards semiflexible mushrooms as the grafting density 
decreases 0- Results from these two coarse-grained 
models are qualitatively the same. Namely, they both 
show the same scaling behavior, but the amplitudes and 
the critical or the crossover points can vary depending on 
the underlying models. However, our main motivation 
was to understand the microscopic differences between 
the two lattice models, SCLM and BFM, for describ¬ 
ing the conformations of polymers and to provide this 
information for the further development of a multi-scale 
approach for studying polymers of complex topology and 
polymer solutions at high concentration based on the lat¬ 
tice models. Therefore, we only focus on the two coarse¬ 
grained lattice models, SCLM and BFM, here. In the 
mapping between atomistic and coarse-grained models, 
it turns out that a bond angle potential also needs to be 
included in the coarse-grained models, as is done here. 

The outline of the paper is as follows: Sec. ITT1 describes 
the models and the simulation technique. Sec. EH and 
Sec. US review the properties of flexible chains and semi¬ 
flexible chains, respectively. Polymer chains described 
by SAWs, and by RWs and NRRWs in the absence of 
excluded volume effects are studied and compared with 
theoretical predictions. Finally our conclusions are sum¬ 
marized in Sec. El 


II. MODELS AND SIMULATION METHODS 

The basic characteristics of linear polymer chains de¬ 
pend on the solvent conditions. Under good solvent con¬ 
ditions the repulsive interactions (the excluded volume 
effect) and entropic effects dominate the conformation, 
and the polymer chain tends to swell to a random coil. 
In the thermodynamic limit, namely the chain length 
N —> oo, the partition function scales as 


(1) 


where /loo is the critical fugacity per monomer, q e g is 
the effective coordination number, and 7 is the entropic 
exponent related to the topology. In two dimensions Q 
7 = 43/32, while the best estimate [El for d = 3 is 
7 = 1.1573(2). For the standard self-avoiding walks on 
the simple cubic lattice [28] in d = 3 one has = 
0.21349098(5) and the corresponding effective coordina¬ 
tion number q e g = l//ioo = 4.6840386(11). The con¬ 
formations of polymer chains characterized by the mean 
square end-to-end distance, (Rj), and the mean square 
gyration radius, (R 2 g ), scale as [29lI30|: 

{Rl)/il = A e N 2v [l + 0{N- A )\, (2) 

and 

{R 2 g)/e l = A g N 2 "{ 1 + 0(N~ a )} (3) 

where v is the Flory exponent, A is the leading correction 
to the scaling exponent, A e and A g are non-universal 
constants, and l 2 is the mean square bond length. The 
quantities u, A, and the ratio A e /A g are universal [3l| . 
while the quantities, A e , A g , t/, and g e ff, depend on the 
microscopic realization. In d = 2 one has 1/2 = 3/4, 
while in d = 3 the most accurate estimate of the Flory 
exponent (32| v = 0.587597(7). We use v = 0.5876 for 
our data analysis in this paper. 

Two models are used for simulating linear polymers in 
the bulk under good solvent conditions. One is the stan¬ 
dard SAW on the simple cubic lattice, effective monomers 
being described by occupied lattice sites, connected by 
bonds of fixed length | b \= £b = 1- Each site can be 
visited only once, and thus the excluded volume inter¬ 
action is realized. The other is the standard bond fluc¬ 
tuation model (BFM). On the simple cubic lattice each 
effective monomer of a SAW chain blocks all eight cor¬ 
ners of an elementary cube of the lattice from further 
occupation. Two successive monomers along a chain are 
connected by a bond vector b which is taken from the 
set {(±2,0,0),(±2,±1,0), (±2,±1,±1), (±2,±2,±1), 
(±3,0,0), (±3,±1,0)} including also all permutations. 
The bond length is therefore in a range between 2 and 
VlO. There are in total 108 bond vectors and 87 dif¬ 
ferent bond angles between two sequential bonds along 
a chain serving as candidates for building the conforma¬ 
tional structure of polymers. The partition sum of a SAW 
of N steps is 


Z N = J2 1 ( 4 ) 

config. 


which is simply the total number of possible configura¬ 
tions consisting of (N ± 1 ) monomers. 

In the literature there are still no estimates of the 
fugacity /i 00 (= 1 /q e fi) and the entropic exponent 7 for 
SAWs on the BFM. According to the scaling law of the 
partition sum Zn, Eq. m, the effective entropic expo¬ 
nent 7 ^ obtained from triple ratios [33| 


Teff (N) = 1 ± 


7 In Zjy — 6 In Zjy /3 — In Z$ n 
ln(3 6 /5) 


Z N ~ ~ g&ATW - 1 


(5) 
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N = L/I b N = L/I b 


FIG. 1: (a) Effective exponents and 7 ® [computed from Eqs. ([5j and (j 6 j ] plotted versus IV on a semi-log scale. 7 = 
limAT-s-oo 7 eff (-2V) = limjv-^oo 7 is (N). = 1.1578(6). (b) In Zn + iVln a — (7 — 1 ) In N with 7 = 1.1578 determined from (a) plotted 
versus N on a semi-log scale. The best estimate of fugacity q, = a = 0.0117241395(75) is determined by the horizontal curve. 


is shown in Fig. QJi. It gives 7 = limjv^oo 7 ^ (TV) = 
1.1578(6). The fugacity \xoa is therefore determined 
by adjusting the value a such that the curve of 
I 11 Zjsr + N In a — ( 7 — l)lnIV with 7 = 1.1578 becomes 
horizontal for very large N (see Fig. We obtain the 
fugacity /ioo = 0.0117241395(75) and the corresponding 
effective coordination number q e g = 85.294106(55) listed 

in Table Q] In Fig. [Q 1 we also show the asymptotic be- 

( 2 ) 

havior of the effective entropic exponent 7/5 defined by 

( 6 ) 

with our estimate of n x for comparison. 

If the excluded volume effect is ignored completely, a 
polymer chain behaves as an ideal chain. It is well de¬ 
scribed by a random walk (RW), a walk that can cross it¬ 
self or may trace back the same path, or by a non-reversal 
random walk (NRRW) where immediate back tracing is 
not allowed. The partition sums of RW and NRRW are 
given by 

Z N = q N (RW), Z N = q(q - l )^” 1 (NRRW) (7) 

where q is the coordination number, q = 6 for the 
standard RW on the simple cubic lattice, and q = 108 
for the BFM [34| . The Flory exponent is v = 1/2 
for an ideal chain and its mean square gyration radius 
(RD = (Rl)/6 . 

For the simulations of single RW, NRRW, and SAW 
chains we use the pruned-enriched Rosenbluth method 
(PERM) [35] . It is a biased chain growth algorithm with 
resampling and population control. In this algorithm 
a polymer chain is built like a random walk by adding 
one monomer at each step with a bias depending on the 
problem at hand, and each configuration carries its own 
weight. The population control at each step is made such 
that the “bad” configurations are pruned with a certain 
probability, and the “good” configurations are enriched 


by properly reweighting, until a chain has either grown 
to the maximum length of steps, N , or has been killed 
due to attrition. A detailed description of the algorithm 
PERM and its applications is given in a review paper |36|. 
The algorithm has the advantage that the partition sum 
can be estimated very precisely and directly. It is also 
very efficient for simulating linear polymer chains up to 
very long chain lengths in dilute solution at and above 
the 0-point. Therefore, we apply the algorithm on the 
two lattice models, SCLM and BFM, in order to check for 
major differences between these two microscopic models. 
The longest chain length is IV = 50000 in our simulations 
here. 


III. CONFORMATIONS OF SINGLE LINEAR 
POLYMER CHAINS: RWS, NRRWS, AND SAWS 

We employ the pruned-enriched Rosenbluth method 
(PERM) for the simulations of long single linear polymer 
chains of chain lengths (segments) up to N ~ 0(1O 4 ), 
modeled by RWs, NRRWs, and SAWs depending on the 
interactions between monomers. Figure [3}i with v = 1/2 
and Fig©) with v = 0.5876 show that the scaling laws, 
Eqs. © and ©, are verified as one should expect. The 
mean square end-to-end distance simply is 



The mean square gyration radius is given by 

"w(S (9) 
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FIG. 2: Mean square end-to-end distance (R 2 ) and gyration radius {Rg) scaled by {N v i 2 ) with the Flory exponent v = 1/2 for 
RWs and NRRWs (a), and v = 0.5876 for SAWs [H}, plotted against N. Here the bond length 4 = (S 2 ) 1 ^ 2 : 4 = 1 (SCLM) 
and 4 = 2.72 (BFM). 
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FIG. 3: Ratio between the mean square end-to-end distance and the mean square gyration radius, ( R 2 )/(R 2 ), plotted against 
the chain segments N, for RWs and NRRWs (a), and for SAWs (b). As N 0, (Re)/(Rg) ~ 6.00 in (a) and {R 2 ) / {R 2 ) ~ 6.25 
in (b). 


where rcM = J2jLo 4/(Af +1) is the center of mass posi¬ 
tion of the polymer. The amplitudes A e and A g for RWs, 
NRRWs, SAWs based on the two lattice models, SCLM 
and BFM, are listed in Table HI Results of {R 2 )/{Nl 2 ) 
and {R 2 )/{Ni 2 ) for RWs from both models follow the 
same curves although the bond vectors in the BFM are 
not all along the lattice directions and do not have the 
same bond length. Here 4 is the root-mean square bond 
length, 4 = 1 for the SCLM and 4 = 2.72 for the BFM. 
In Fig. [2}a, values of {R 2 )/{N£ 2 ) and {Rg) /{Ni 2 ) for NR¬ 
RWs, obtained from SCLM for all lengths N are sig¬ 
nificant larger than that obtained from the BFM, since 
at each step the walker can only go straight or make 
a 90° L-turn in the SCLM. In Fig. [2]o, the two curves 
showing the results of {R 2 e ) / (N v t 2 b ) [(R 2 g ) / (N v £ 2 )] with 
v = 0.5876 as functions of N obtained from the two mod¬ 
els intersect at N « 180, and finally the amplitude for 
BFM is larger in the asymptotic regime. The slight de¬ 
viation from the plateau value is due to the finite size 


effects. The correction exponent A [Eqs. © and ©] for 
these two models is determined by plotting {R 2 )/N 2v and 
{R 2 )/N 2v versus x = N~ s (not shown). One should ex¬ 
pect straight lines near x = 0 if and only if <5 = A. We ob¬ 
tain A = 0.48(5) for both models, which is in agreement 
with the previous simulation in Ref. [29|, [32j, [H] within er¬ 
ror bars. The ratio between the mean square end-to-end 
distance and the mean square gyration radius, is indeed 
{R 2 e )/{R 2 g ) « 6.000(3) for RWs and NRRWs (Fig. Et)- 
For SAWs our results give {R 2 e )/(R 2 g ) = 6.25(2) (Fig.EJa). 
For SAWs on the simple cubic lattice the most accurate 
estimates of A e = 1.22035(25), A g = 0.19514(4), and 
A e /A g = 6.2537(26) are given in Ref. [32j. Our results 
are also in perfect agreement with them. However, much 
longer chain lengths will be needed for a more precise es¬ 
timate of the plateau value of the ratio in the asymptotic 
scaling regime. Note that for 10 < N < 100 the behavior 
is clearly model-dependent. 

We include here the RW and NRRW versions of both 
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TABLE I: The estimates of fugacity /ioo, the effective coordination number g e ff in Eq. ©, the amplitudes A e and A g in Eqs. © 
and m determined from the simulation data of (R 2 ) and (R 2 g ) for RWs, NRRWs, and SAWs based on the two lattice models, 
SCLM and BFM. 



SCLM 

BFM 

model 

RW 

NRRW 

SAW 

RW 

NRRW 

SAW 

fJ'OO 

1/6 

1/5 

0.21349098(5) [28] 

1/108 

1/107 

0.01172414395(75) 

QeS 

6 

5 

4.6840386(11) 

108 

107 

85.294106(55) 

A e 

1 .0000(2) 

1.4988(4) 

1.220(3) 

0.9986(2) 

1.0714(2) 

1.247(5) 

A g 

0.16666(0) 0.24985(7) 

0.1952(4) 

0.16645(4) 0.16959(3) 

0.1993(6) 
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FIG. 4: (a),(c) Normalized probability distributions of end-to-end distance, hN(Re/£b)(= V(Re/ib)), plotted versus R e /£b- 
(b),(d) similar as (a),(c), but for gyration radius R g . Data are for RWs (a),(b) and SAWs (c),(d). Several values of chain 
lengths N are chosen, as indicated. 


models not just for the sake of an exercise: often the 
mapping from an atomistic to a coarse-grained model is 
to be done under melt conditions, where excluded volume 
interactions are screened. 


The shapes of polymer chains can also be described by 
the probability distributions of end-to-end distance and 
gyration radius, P(R e /£b) and P(R g /£b), respectively. 
Numerically, they are obtained by accumulating the his¬ 
togram Hn(x) of x over all configurations of length N, 


given by 

H n (x)= Y, W N (x')S x , x ', (10) 

configs. 

here each configuration carries its own weight Wn(x'). 
The normalized histogram is therefore, 

tiN(x) = H n (x)/ Y H n (x') . (11) 

x' 

Results of h N (R e /i b ) and h,N(R g /h) for RWs and SAWs 
obtained from the two models for various values of chain 
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FIG. 5: Same as in Fig. [dja,, but data are only for BFM. The 
predicted distribution 4'K(R e /£b) 2 P(R e /£b) for various values 
of N are shown by solid curves. Here the distribution function 
P{R e /£b) is given by Eq. (fl6l) . 



FIG. 6: Root mean square radius of gyration, (R 2 ) 1 ^ 2 , and the 
gyration radius R g ,m where P(R g ) has its maximum value, 
plotted against chain length N for RWs and SAWs. Data are 
for BFM. 


lengths N are shown in Fig. [4] We see that both mod¬ 
els give for N > 100 the same distributions of R e /h 
and Rg/ib although the mean values of R e /h and R g /£ b 
are slightly different between these two lattice models 
(Fig. [Do) for SAWs. Note that an angular average over 
all directions has been included in the accumulating pro¬ 
cess of the histogram due to spherical symmetry. Thus, 
the normalized histograms of R e /£b , 

h N {Re/h) = V N {Re/ib) = 4 nC e , N (-R e /4) 2 P N ($*/h) , 

( 12 ) 

with 

T N (R e /£ b )d(R e /£ b ) = l, (13) 

and the normalized histograms of R g /£ b , 

h N ( R g /h ) = V N {Rg/£b) = 4 TrC g , N ( Rg/ibf P N (R g / 4) , 

(14) 

with 

V N (Rg/£ b )d{Rg/e b ) = l, (15) 




conditions in d-dimensions suggested by Lhuillier [37] is 
as follows: 


P(R g /£ b ) ~ exp 


-a i 


(£ b N v \ ad ( R g \ 
\ Pg ) ° 2 {hN") 


(17) 

where ai and 02 are (non-universal) constants, and the 
exponents a and 5 are linked to the space dimension d 
and the Flory exponent v by 


a = ( vd — 1) 1 and S = (1 — v) 1 . (18) 

Here (1 + a) is the des Cloizeaux exponent (38] for the 
osmotic pressure of a semidilute polymer solution, and <5 
is the Fisher exponent [39j characterizing the end-to-end 
distance distribution. 

This scaling form has been verified in the previous 
Monte Carlo simulation studies of the standard self¬ 
avoiding walks on square ( d = 2) and cubic (d = 3) 
lattices up to O(l0 2 ) step s using the slithering-snake and 
pivot algorithms [40, [4l|. The two fitting parameters ai 
and 02 are actually not independent since at the position 
where the distribution P(R g ) has its maximum value, i.e. 
P(Pg = Pg,m) = max P(R g ), the corresponding gyra¬ 
tion radius -R s , m oc R g oc N" (see Fig. [UI)- Using Eq. (fT7[i . 
the logarithm of the rescaled probability is written as 


where C^jv and C g are the normalization factors. 

The probability distribution of end-to-end distance for 
ideal chains is simply a Gaussian distribution, 

(i 6 ) 

Our numerical data for RWs obtained from BFM and 
SCLM shown in Fig. [4] are in perfect agreement with the 
Gaussian distribution (see Fig. 0 . From Eqs. (usd, (usd 
and (1131) we obtain the normalized factor C e ,n = 1. 

The theoretical prediction of the gyration radius prob¬ 
ability distribution of polymer chains under good solvent 


with 


a 1 = 


f 

= A 

A f R, 


Ro 


Rg 

1 (R a 


= In 


P{Rg,m/h) 
P(Rg/h ) 


R„ 


ad 


Ra 


R 


■g,m 


1 — d 


a \e b N" 

From Eq. asD, we obtain 


and 02 = 


Ad (£ b r 


R 


■g,m 


P(Rg,m/£b) _ hN(Rg,m/£b) _|_ 2 ^ 


P(Rg/h ) 


hN(Rg/h ) 


R, 


g,m 


(19) 


( 20 ) 


( 21 ) 
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FIG. 7: Logarithm of the rescaled probability distribution of gyration radius, In {P(R g ,m) / P{R g )) as a function of R g /R g ,m. 
for SAWs obtained from the models (a) SCLM and (b) BFM. The best fit of Eq. (11911 with A = 1.18 to our data is shown by 
the solid curve. The dashed curve with A = 1.34 given in Ref. [TL] is also shown for the comparison. Several values of chain 
lengths N are chosen, as indicated. 
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FIG. 8: Similar as in Fig. [7] but for RWs. The best fitting of our data gives A = 0.97 for both models. 


Our estimate of ln(P(R gtrn /£b)/P(Rg/h)) for SAWs 
based on the two lattice models, SCLM and BFM, are 
shown in Fig. [7] As chain lengths N > 1000, we see the 
nice data collapse, and the logarithm of the scaled prob¬ 
ability of R g is described by Eq. © with A = 1.18 very 
well. Due to the finite-size effect it is clearly seen that 
the previous estimate A = 1.34 is an overestimate [4lj . 

For an ideal chain the distribution of R g is no longer 
a simple Gaussian distribution as shown in Eq. ©, and 
the exact expression is quite complicated. Vettorel et 
al. 0 found out the formula given by Lhuillier [37]] 
is a good approximation for describing the distribution 
P{R g ) for an ideal chain based on the BSM. Therefore, 
we also use the same formula for the investigation of the 
distribution P(R g ) obtained from the two coarse-grained 
lattice models. Two methods are discussed here. Method 
(1): We use the formula [Eq. (fl9l) ] which contains only 
one fitting parameter A since R g , m ~ R g ~ N for RWs 
as shown in Fig. [G] From our simulations of RWs, we 
still see the nice data collapse for N > 200 in the plot of 


the logarithm of the rescaled distribution of R g (Fig. [S]), 
but the distribution can only be described by Eq. © 
well for R g > R gjrn . Using the least square fit, it gives 
A = 0.97. Method(2): We assume that the two pa¬ 
rameters a\ and <22 in Eq. © are independent. Us¬ 
ing Eqs. ©, ©, and ©. values of ai, < 22 , and the 
normalization factor C 9t jv are determined by the best fit 
of the normalized histograms h]y(R g /£b ) obtained from 
our Monte Carlo simulations. Note that it is not pos¬ 
sible to determine a\ and 02 using the second method 
for N < 50 since the normalization condition, Eq. ©, 
is not satisfied. In Fig. [9] we compare our results of 
h N (R g /£b ) OC V{Rg/h) OC ( R g /h) 2 P{Rg/e b ) for BFM to 
the fitting function 47 t C g N{R g /£b) 2 P{R g /h) [Eqs. (070) . 
©, and (TT7T) ] with parameters determined by these two 
different methods. Values of a\ and a 2 plotted versus N 
are shown in Fig. [TO] and listed in Table El Our results 
show that ai and 02 are almost constants for large N, 
which are comparable with the results obtained for the 
BSM 0. 























(a) 


(b) 




FIG. 9: Normalized probability distributions of R g /£b, h{R g /lb) = VN(R g /£b), plotted versus R g /lb for various values of A, 
and for BFM. The fitting functions 47rCg,jv (. R a /£ b ) 2 P{R g /£b ) [Eqs. (HI, US) , and ED] with parameters £2i, 02 ■ and C gi N 
determined by method (1) and method (2) are shown by curves for comparison in (a) and (b), respectively. 


TABLE II: Parameters ai and a 2 of the probability distribution P(R g /£b), Eq. (1171) . determined by two different methods (1) 
and (2) mentioned in the text for various values of chain length N. 


N 

10 

20 

50 

100 

200 

400 

800 

1000 

2000 

4000 

8000 

10000 

(1) ai x 10 4 

7.12 

5.27 

4.94 

4.52 

4.70 

4.70 

4.82 

4.34 

4.88 

4.74 

4.37 

4.68 

(!) °2 

12.80 

14.15 

14.46 

14.90 

14.70 

14.73 

14.58 

15.10 

14.42 

14.66 

15.07 

14.73 

(2) ai x 10 4 

_ 

_ 

4.14 

3.83 

3.69 

3.60 

3.57 

3.56 

3.58 

3.54 

3.53 

3.53 

(2) <22 

_ 

_ 

13.25 

13.27 

13.29 

13.29 

13.30 

13.29 

13.35 

13.29 

13.29 

13.28 


1.6 
1.5 
1.4 
1.3 
1.2 
1.1 
1 

0.9 
0.8 

10 100 1000 10000 
N 

FIG. 10: Parameters a\ and <22 in Eq. (H3 plotted versus N. 
Results are obtained from two different methods mentioned 
in the text. Here ai i00 and a 2 , 0 a are taken from Table HT1 for 
N = 10000. 


IV. SEMIFLEXIBLE CHAINS 

We extend our simulations in this section from flexi¬ 
ble chains to semiflexible chains. Extensive Monte Carlo 
simulations of semiflexible polymer chains described by 
standard SAWs on the simple cubic lattice, with a bend¬ 
ing potential U b = e b (l — cos 9), have been recently car¬ 


ried out [42l - l44l |. Recall that atomistic models of real 
chains may exhibit considerable chain stiffness due to 
the combined action of torsional and bond angle poten¬ 
tials. When a mapping to a coarse-grained model is per¬ 
formed, this stiffness is lumped into an effective bond 
angle potential of the coarse-grained model. In this stan¬ 
dard model the angle between two subsequent bond vec¬ 
tors along the chain is either 0° or ±90°, and hence in 
the statistical weight of a SAW configuration on the lat¬ 
tice every 90° bend will contribute a Boltzmann factor 
qb = exp(— £b/ksT) (q b = 1 for ordinary SAWs). fcsT is 
of order unity throughout the whole paper. The partition 
function of such a standard SAW with A bonds (A + 1 
effective monomers) and Abend bends is therefore, 

Z N (q b )= ^v(Abend)gf bend , (22) 

config. 

where Cat (Abend) is the total number of all configurations 
of a polymer chain of length A containing Abend bends. 

We are also interested in understanding the micro¬ 
scopic difference between the standard SAWs and the 
BFM as the stiffness of chains is taken into account. Since 
there are 87 bond angles possibly occurring in the chain 
conformations, the partition function cannot be simpli- 


a i / a i,« 
a i / a i 

/ &2 oc 
/ &2 oc 

method (2) — b — 
method (1) 
method (2) * 

method (1) 

k. \ 


"X 


- -Sr-"' * 
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fied for the BFM, written as, 


4 bfm) m 

= Cat({6 »}) exp 

config. 


N-l 
£b ST' 


knT ^ 


(1 - cos Oi,i+i) (2,3) 




where Oi^+i is the bond angle between the i th bond vec¬ 
tor and the (i + l) th bond vector along a chain, and 
Cat({0 }) is the number of configurations having the same 
set {0} but fluctuating bond lengths. In the absence 
of excluded volume effect, the formulas of the partition 
function, Eq. (T221) and Eq. (E51) . remain the same while 
semiflexible chains are described by RWs and NRRWs. 


A. Theoretical predictions 


There exist several theoretical models describing the 
behavior of semiflexible chains in the absence of excluded 
volume effects. We first consider a discrete worm-like 
chain model [45[ that a chain consisting of N bonds with 
fixed bond length £b, but successive bonds are correlated 
with respect to their relative orientations, 

(bi ■ b i+ 1 ) = ^(cos 9) and (b 2 ) = £ 2 h , (24) 

where 9 is the angle between the successive bond vectors. 
The mean square end-to-end distance is therefore, 

l + (cos0) 2(cos0)(l — ((cos9)) N 
1 — (cos 9) N( 1 —(cos0)) 2 

(25) 

This formula agrees with the prediction for a freely ro¬ 
tating chain (FRC). In the limit N —> oo the bond-bond 
orientational correlation function decays exponentially as 
a function of their chemical distance s, 

(bi ■ b i+s ) = £l(cos9(s)) = (l(cos9) s = £ 2 exp (-s £ b /£ p ), 

(26) 

where £ p is the so-called persistence length which can be 
extracted from the initial decay of (cos 9(s)). Equiva¬ 
lently, one can calculate the persistence length from 


(Rl) = m 2 b 


£ P ,s/£b = -1/ ln((cos 9)) (27) 


here instead of £ p we use £ Pi $ to distinguish between these 
two measurements. 

For rather stiff (L £ P ) and long chains (N —> oo) we 
expect that the bond angles 9 between successive bonds 
along chains are very small ((cos 9) « 1 — (9 2 )/2 ), then 
Eqs. m and m become 


(R 2 e ) = m 2 b 


1 + (cos 9) 
1 — (cos 9) 



(28) 


and 


Eq. (1281) is equivalent to the mean square end-to-end dis¬ 
tance of a freely jointed chain that rip Kuhn segments of 
length ix are jointed together, 

(Rl)=n k tl = 2i p L. (30) 

L = Ntp = n k £k being the contour length and £k = 2£ p 
in this limit. 

In the continuum limit £b — > 0, N —> 0, but keeping L 
and £ p finite, we obtain from Eq. (1251) the prediction for 
a continuous worm-like chain, 


(Rl) = 2 £ P L { 1 - j-[l - exp (~L/ip)\ 


(31) 


It gives the same result as that derived directly from 
the Kratky-Porod model [H, j4l} for worm-like chains in 


d = 3, 


t P k B T f L / d 2 r(s) 

2 J 0 1 ds 2 


2 

ds , 


(32) 


where the polymer chain is described by the contour 
r(s) in continuous space. Equation (fSTI) describes the 
crossover behavior from a rigid-rod for L < £ p , where 
(Rl) = L 2 , to a Gaussian coil for L £ p , where 
(Rl) = 2 £ p L as shown in Eq. (1301) . 

For semiflexible Gaussian chains the contour length 
L = N£b can also be written as L = n p £ p and the mean 
square end-to-end distance and gyration radius described 
in terms of n p and £ p are [46i . |48[ 


{Rl) 

2 £ p L 


= 1-[1 — exp(— n p 


(33) 


6 (R 2 g ) _ x 3 | 6 
2£ P L Tip 


—exp(-n p )]. 
n p 


(34) 


One can clearly recognize that Gaussian behavior of the 
radii is only seen, if the number n p of the persistence 
length that fits to a given contour length is large, n p 1, 
while a crossover to rigid-rod behavior occurs for n p of 
order unity. 

In recent works in Ref. J2^, 43, 43| , authors have shown 
that the exponential decay of the bond-bond orienta¬ 
tional correlation function, Eq. (l26l) . and the Gaussian 
coil behavior, Eq. m for L £ p , predicted by the 
wornr-like chain model only hold for s and N up to 
some values s* and N *, respectively when excluded vol¬ 
ume effects are considered. The predictions of a theory 
based on the Flory-type free energy minimization argu¬ 
ments @, MM proposed as an alternative to semiflex¬ 
ible chains with excluded volume interactions have been 
verified. In this treatment one considers a model where 
rods of length £ k and diameter D are jointed together, 
such that the contour length L = N£b = n k £k■ Apart 
from prefactors of order unity, the second virial coeffi¬ 
cient in d = 3 then can be estimated as 


(29) 


= £\D. 


£ P /£ b = 2 /(9 2 ) 


(35) 
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The free energy of a chain now contains two terms, the 
elastic energy taken as that of a free Gaussian, i.e., F e i ss 
R 2 /{£ k L), and the repulsive energy due to interactions 
treated in mean field approximation, i.e. proportional 
to the square of the density n/R 3 and the volume R 3 . 
Hence, 

A F/k B T « R 2 J(£ k L) + v 2 Rl[(L/i K )/Rlf (36) 

Minimizing A F with respect to R e , we obtain for L —> oo 
the standard Flory result 

Re » («i/4) 1 / 5 L 3/5 = (4£)) 1 / 5 (AT4 ) 3/5 • (37) 

Eq. (l37l) holds also for finite L and N > N* = £\/{£bD 2 ) 
since the contribution of the second term in Eq. (1361) 
is still important. For N < N* the first term in 
Eq. (1351) dominates, and the chain behaves as a Gaus¬ 
sian coil, R 2 = £ k L = £ k £bN, while for even smaller N , 
N < N rod = 4/4, the chain behaves as a rigid-rod. 
Thus, the double crossover behavior of the mean square 
end-to-end distance is summarized as follows, 

(R 2 ) ss L 2 , N < N md = 4/4 (rod — like chain), 

(38) 

(R 2 e ) « 4 L , N rod < N < N* (Gaussian coil), (39) 
(R 2 e ) ss ( 4 D) 2 / 5 L 6 / 5 , N > N* (SAW) (40) 


B. Simulation Results 

In order to investigate the scaling behavior of the ratio 
(R 2 )/(R 2 ) for semiflexible RWs and NRRWs we plot our 
data (R 2 )/ {R 2 g ) versus N = L/4 for several choices of the 
stiffness parameter (Fig. 1TT1) . As N increases, the data 
increase towards a maximum and then decrease towards 
a plateau where the prediction limjv-^oo {R 2 )/(R 2 g ) ~ 6 
for ideal chains holds. At the location of the maxi¬ 
mum of ( R 2 )/{R 2 ), N = Nh, the corresponding maxi¬ 
mum values is h = max (R 2 )/(R 2 ). The maximum move 
monotonously to larger values as chains become stiffer. 
The deviation between the data for RWs and NRRWs 
based on the SCLM decreases as the bending energy Eb 
increases (Fig. [lTJi), while it is negligible for the simula¬ 
tion data obtained based on the BFM (Fig. ITTbl in all 
cases. 

Figure [T5] shows the bond-bond orientational corre¬ 
lation function (cos (9(s)} plotted versus the chemical 
distance s 4 covering the range from flexible chains to 
stiff chains characterized by Eb for the models SCLM 
and BFM. We compare the data obtained for SAWs, 
NRRWs, and RWs for various values of Eb- The in¬ 
trinsic stiffness remains the same for SAWs, NRRWs, 
and RWs as Eb is fixed. Results obtained from both 
models verify that the asymptotic exponential decay of 


(cos 9(s)) is valid only if the excluded volume effect is 
neglected, i.e., for RWs and NRRWs. For semiflexible 
SAWs (cos0(s)) ~ exp(—s4/4) cannot be correct for 
N —> oo [22|, [42J , we rather have 

(bi ■ b i+s ) « s-P , P = 2 — 2v S3 0.824 , s* < s < N . 

(41) 

As we have seen in Fig. [ 12 l the exponential decay is 
ill-defined for rather flexible SAWs. Using Eq. (1571) 
as a definition of the persistence length we can still 
give the estimate of the persistence length £ Pt g = 
— 4 /In [(cos 0 (s = 1 ))] which is approximately the same 
as the estimate of the decay length £ p for moderately stiff 
chains and stiff chains. The estimates of £ p /4 and £ P: g/£b 
depending on Eb using Eqs. (051) and (071) are listed in Ta¬ 
ble [TTIl and IIVI RWs are more flexible than NRRWs, and 
NRRWs are more flexible than SAWs from the estimates 
of the persistence lengths £ p /4 and £ Pt g/£b based on the 
SCLM. Using the BFM the persistence lengths are almost 
the same in all cases of Eb for RWs and NRRWs, and 
they are smaller compared with the estimates for SAWs. 
Note that in Fig. [T2b data deviate slightly from the fit¬ 
ting straight lines describing the initial exponential decay 
for RWs and NRRWs as the bending energy Eb increases, 
i.e., the stiffness of chains increases. For Eb > 10 the 
problem is more severe. Therefore, one should be careful 
using the BFM for studying rather stiff chains. An alter¬ 
native way to the determination of the persistence length 
would be given by the best fit of the mean square end- 
to-end distance {R 2 ) of RWs or NRRWs to Eq. (1311) . A 
simple exponential decay is always found for the prob¬ 
ability distribution of connected straight segments for 
semiflexible chains based on the SCLM [43j], while large 
fluctuations are observed for semiflexible chains based 
on the BFM due to bond vector fluctuations and lattice 
artifacts [52|. This is the main reason why the different 
scenarios of the bond-bond orientational correlation func¬ 
tions between the SCLM and the BFM for stiff chains are 
seen in Fig. [T7] Figure [13] shows the locations Nh and 
the heights h of the maximum of {R 2 )/{R 2 ) (Fig. fill) 
plotted versus the persistence length £ p /4 for semiflexi¬ 
ble RWs based on the SCLM and BFM. Note that Nh, h , 
and £ p /£b all depend on £5 which controls the stiffness of 
chains. We see that the dependence between h and £ p /4 
are the same for both models, while N k for the BFM is 
slightly larger than that for the SCLM for a fixed value 
of £ p /£b since chains based on the BFM are more flexible. 

The scaling plots for testing the applicability of the 
worm-like chain prediction, Eq. (IBTl) and Eq. (1551) to our 
data of (R 2 ) and (R 2 ) are shown in Fig. [15] The persis¬ 
tence length 4/4 in Eq. (1511) for various values of Eb are 
extracted from the exponential fit of Eq. (1551) for NRRWs 
(see Tables IIIII and IIVI) . Since the wornr-like chain model 
is formulated in the continuum, care has to be taken to 
correctly take into account the lattice structure of the 
present model, particularly in the rod limit. Assuming 
that a rigid rod consisting of N monomers is located at 
x\ = 4, £2 = 24, • • Xn = N along the x-axis on the 
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(a) (b) 




FIG. 11: Ratio between mean square end-to-end distance and mean square gyration radius, (R 2 )/(Rg), plotted against chain 
lengths (segments) N = L/lb for SCLM (a) and for BFM (b). Data for semiflexible chains described by NRRWs and RWs are 
shown by symbols and lines, respectively. Here ih = 1 for SCLM, and lb = 2.72 for BFM. 


(a) (b) 




FIG. 12: Semi-log plot of the bond-bond orientational correlation function (cos 9(s)) vs. sib for SCLM with lb = 1 (a) and for 
BFM with lb = 2.72 (b). Data are for semiflexible chains described by SAWs, NRRWs and RWs and for e b = 5.30, 4.61, 3.91, 
3.51, 3.00, 2.30 from top to bottom in (a), and for Eb = 10, 7, 5, 3, and 1 from top to bottom in (b). The straight lines indicate 
fits of the initial decay, (cos#(s)) oc exp(— slb/l P ) [Eq. (1261) ]. for RWs. 


simple cubic lattice, the mean square gyration radius is: 

I N ( i N \ 2 

(-Re)rod = — y^(fclfc ) 2 - I TF fclfr ) 
k= 1 \ k= 1 / 

(N+1)(2N + 1)£ 2 (N+l) 2 £ 2 

6 4 

(N + 1)(N-1)£ 2 _ L(L + 24) , . 

12 12 ' 1 ’ 

Therefore, due to the lattice structure, the mean square 
gyration radius is rescaled by (L + 2 lb) instead of L 
in order to compare with the theoretical predictions in 
Fig. H4b.d. For semiflexible RWs and NRRWs the data 
are indeed very well described by the worm-like chain 
model. As N increases, we observe the crossover behav¬ 
ior from a rigid-rod regime to a Gaussian coil regime. 
The plateau value in the Gaussian regime corresponds to 
the persistence length £ p /£b in Fig. [LUt.b and (1/6 )£ p /£b 


in Fig. ITdb.d. For SAWs the deviation from the predic¬ 
tion becomes more prominent as chains are more flexible 
since the excluded volume effects are more important. 

Note that one should not consider the correction fac¬ 
tor (L + 2 £b)/L relative to the Kratky-Porod model in 
Eq. (H51) as a “lattice artefact”: in a real stiff polymer 
(e.g. an alkane-type chain) one also has a sequence of 
discrete individual monomers (separated by almost rigid 
covalent bonds along the backbone of the chain) lined 
up linearly (like in a rigid rod-like molecule) over about 
the distance of a persistence length. Furthermore, we 
compare simulation results of the ratio (R 2 )/(R 2 ) mul¬ 
tiplied by \(L + 2 £b)/L\ as a function of n p = L/£ p to 
the theoretical prediction, the ratio between Eq. (1531) 
and Eq. (l34l) . in Fig. [TS] We see the nice data collapse 
for RWs and NRRWs in the Gaussian regime ( n p 1) 
and the increase of deviations from the master curve as 
the stiffness of chains decreases in Fig. flSli.b. The ra- 
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TABLE III: Two estimates for the persistence length £ p /£b from Eq. (1261) and 4,9/4 from Eq. (1271) for semiflexible RWs, 
NRRWs, and SAWs with various values of qb{= exp(—et/fesT)) based on the SCLM ( £b = 1). Here in our simulations values 
of qb are chosen for convenience. 



96 

1.0 

0.4 

0.2 

0.1 

0.05 

0.03 

0.02 

0.01 

0.005 


Eb 

0.0 

0.91 

1.61 

2.30 

3.00 

3.51 

3.91 

4.61 

5.30 

4/4 

RW 


0.84 

1.54 

2.83 

5.37 

8.73 

12.95 

25.67 

51.38 


NRRW 


1.05 

1.70 

2.97 

5.50 

8.87 

13.09 

25.80 

51.53 


SAW 



2.04 

3.35 

5.96 

9.54 

13.93 

26.87 

52.61 

4 . 9/4 

RW 


0.84 

1.54 

2.83 

5.36 

8.73 

12.95 

25.66 

51.37 


NRRW 

0.62 

1.05 

1.70 

2.98 

5.50 

8.87 

13.08 

25.79 

51.53 


SAW 

0.67 

1.12 

1.81 

3.12 

5.70 

9.10 

13.35 

26.28 

51.52 


TABLE IV: Two estimates for the persistence length £ p /£b from Eq. (1261) and 4,9/4 from Eq. (1271) for semiflexible RWs, 
NRRWs, and SAWs with various values of £b based on the BFM (£b = 2.72). 



Eb 

0.0 

1.0 

2.0 

3.0 

5.0 

7.0 

10.0 

15 

4/4 

RW 


0.87 

1.62 

2.54 

4.69 

7.18 

12.09 

27.65 


NRRW 


0.87 

1.62 

2.54 

4.69 

7.18 

12.09 

27.65 


SAW 



1.91 

2.78 

4.94 

7.39 

12.37 

27.93 

4 . 9/4 

RW 


0.87 

1.62 

2.54 

4.63 

6.87 

10.50 

17.73 


NRRW 

0.21 

0.87 

1.62 

2.54 

4.63 

6.87 

10.50 

17.73 


SAW 

0.61 

1.11 

1.80 

2.65 

4.68 

6.90 

10.52 

17.75 



FIG. 13: Location Nh and height h of the maximum of 
(Re)/{Rg ) (see Fig. 1111) plotted versus the persistence length 
4/4 for RWs based on the SCLM and BFM. 


tio [(L + 24 )/L](-Rg)/(i?g) ~ 12 as n p —> 0 for a rigid- 
rod, while [(L + 2£b)/L\(Rl)/(R^) ss 6 as n p —> oo for a 
Gaussian coil. For SAWs we still see the nice data col¬ 
lapse in Fig. [TSb.d. but in both rigid-rod and Gaussian 
coil regimes the deviations from the master curve become 
more prominent as chains are more flexible. For n p > 1 
the deviation is due to the excluded volume effects, and 
finally [(L + 2£j ,)/L\(Rl)/(R^) « 6.25 as n p —> oo for 
SAWs. Note that in both models the ratio of the mean 


square end-to-end and gyration radii exceed its asymp¬ 
totic value still significantly even if n p is as large as 
n p = 10. 

Recently, Huang et al. [H],[54| performed Brownian dy¬ 
namics simulations on two-dimensional (2D) semiflexible 
chains described by a BSM including the excluded vol¬ 
ume interactions. Varying the chain stiffness and chain 
length their results confirmed the absence of a Gaussian 
regime in agreement with the results from semiflexible 
SAWs based on the SCLM jd4|, and with observations 
from experiments of circular single stranded DNA ad¬ 
sorbed on a modified graphite surface [? ]. The rescaled 
mean square end-to-end distance, (Rl)/(2L£ p ), in terms 
of L/i v for both models on the lattice and in the contin¬ 
uum turns out to be universal from the rigid-rod regime 
up to the crossover regime ( L/£ p ~ 1) irrespective of 
the models chosen for the simulations. In the 2D SAW 
regime, different amplitude factors result from the differ¬ 
ent models [25| . 

In d = 3, we indeed see the nice data collapse for 
semiflexible RWs, NRRWs, and SAWs in the plot of 
(Rl)/(2L£ p ) versus L/2£ p (cf. Fig. fLlh.b) from rod-like 
regime crossover to the Gaussian regime for N < N* 
(not shown), and the data obtained from the two lattice 
models are well described by the Kratky-Porod scaling 
function, Eq. m- For the BSM in the continuum we 
should expect the same universal behavior. Although for 
semiflexible SAWs the second crossover from the Gaus- 
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(a) 


(b) 


A 

CM d) 

DC 
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(c) (d) 




FIG. 14: Log-log plots of rescaled mean square end-to-end distance (R 2 )/(2£bL) (a),(b) and rescaled mean square gyration 
radius (R 2 ) / (2£b(L + 2£b)) (c),(d) versus N = L/£b for semiflexible chains described by SAWs, NRRWs and RWs based on the 
SCLM with £b = 1 (a),(c) and BFM with £b = 2.72 (b),(d). Data for various values of Eb are shown, as indicated. Solid curves 
refer to the theoretical prediction, Eq. (O, for WLC. The values of the persistence length £ p /£b for NRRW are taken from 
Tables Hill and Hvl 


sian regime to the SAW regime for N > N* is rather 
gradual and not sharp, the relationship [5fil | between the 
crossover chain length N* and the persistence length 
£ p /£b , A* oc (ip/h) 2 ' 5 , holds for these two models here. 
It would be interesting to check whether such a scaling 
law would also hold for the BSM. 

The structure factor S(q) is an experimentally accessi¬ 
ble quantity measured by neutron scattering. We there¬ 
fore also estimate S(q) by 


S(q) 


1 


(JV + l ) 2 


! N N 

EE exp(*<f • \ri 

\ 2—0 j — 0 



(43) 


where (f)} denote the positions of the (iV + 1) monomers 
in a chain, and the structure factor is normalized such 
that S(q —> 0) = 1. In order to compare the results of 
S(q) obtained for fully flexible RWs, NRRWs, and SAWs 
based on the SCLM and BFM, we plot S(q) versus q£b 
(4 = 1 for SCLM, and i b = 2.72 for BFM) in Fig. [[Sl. 
We see that S(q) « 1 — q 2 (R 2 )/3 for q —> 0, while for 

q y/ (Rg) the power law S(q) ~ q~ x ! v (v = 0.588 for 

SAWs, and v = 0.5 for RWs and NRRWs) holds. The lat¬ 
tice artifact sets in at qi b ~ i r. Due to the local packing 


the first peak appears at qi b ~ 2-7T for the SCLM, while 
at q£b ~ 2.47T for the BFM as q increases. In Fig. 17771 ) we 
show the results for semiflexible SAWs of different stiff¬ 
nesses based on the BFM. The Gaussian regime where 
S(q) ~ q~ 2 for large values of £ b and then crosses grad¬ 
ually over to S(q) ~ q~ x as expected for rigid rods. [57j . 

Finally we analyze the structure factor S{q) in the 
form of Kratky-plots, qLS(q) plotted versus ql p , shown in 
Fig. 1777 for semiflexible SAWs. Data are only for q < n. 
The well-known theoretical predictions of the scattering 
from rigid-rods [§3], qLS(q) —>• 7 r, and Gaussian chains, 
the Debye function |j. |58h6G ] , 

= T^tf) {'- iFS !) 11 -, 

(44) 

and the interpolation formula which describes the two 
limiting cases of Gaussian coils and rigid rods exactly by 
Kholodenko [ 6 l| . 

S{q) = -\h(x) - -J 2 (x)], x = 3L/2l p (45) 
xx 

x 

where I n (x) = / dzz 71-1 f(z), and the function f(z) is 
o 
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(a) (b) 






FIG. 15: Semi-log plots of [(L + 2 £b) / L\{R 2 ) / (R 2 ) versus n p = L/£ p for semiflexible chains described by RWs and NRRWs 
(a),(b) and SAWs (c),(d). Data for various values of Eh are shown, as indicated. Solid curves refer to the theoretical prediction, 
the ratio between Eq. (1331) and Eq. (1341) . for WLC. The values of the persistence length £ p /£b for RWs, NRRWs, and SAWs are 
taken from Tables m and IIVI respectively. 


given by 


/(*) 


1 sinh(gz) < _3_ 

E sinh z ** — 2 t v ’ 
1 sin(gg) ^ 3 

sinh 2 5 ** 2 i v 7 


(46) 


with 



(47) 

are also shown for comparison [62|. Near the peak the 
discrepancy of our data from the theoretically predicted 
formulas increases as the bending energy E\, decreases 
showing that the excluded volume effect sets in. For 
semiflexible polymer chains of almost the same persis¬ 
tence length based on the two different lattice models, 
the structure factors are on top of each other. 


V. CONCLUSIONS 

In this paper we have studied single polymer chains 
covering the range from fully flexible chains to stiff chains 


under very good solvent conditions by extensive Monte 
Carlo simulations based on two coarse-grained lattice 
models: the standard simple cubic lattice model and 
the bond fluctuation model. With the pruned-enriched 
Rosenbluth method the conformations of polymer chains 
mimicked by random walks, non-reversal random walks, 
and self-avoiding walks depending on the effective inter¬ 
actions between monomers have been analyzed in detail. 
We give the precise estimate of the fugacity [Xoo and the 
entropic exponent 7 for self-avoiding walks based on the 
bond fluctuation model. The universal scaling predic¬ 
tions of mean square end-to-end distance, {R 2 e ) [Eq. ©], 
and mean square gyration radius, (R g ) [Eq. Q], for fully 
flexible chains are verified as one should expect, and the 
corresponding amplitudes A e and A g depending on the 
models are determined. We have also checked the prob¬ 
ability distributions of R e and R g , P(R e ) and P(R g ) 7 
respectively. Especially we point out that the previous 
estimate of the parameter A in Eqs. (HU for SAWs is 
an overestimate due to the finite-size effect. Our results 
also agree with the results based on the BSM 0 , that 
the formula Eq. (ED predicted by Lhuillier [371) is a good 
approximate formula for RWs. 

For semiflexible chains the additional regime of rod-like 
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(a) (b) 




q'b 


q'b 


(c) 



FIG. 16: (a)(b) Log-log plot of structure factor S(q) versus qib . Data are for fully flexible chains based on the SCLM and 
BFM in (a), and for semiflexible chains based on the BFM including 6 choices of the stiffness in (b). (c) Rescaled structure 
factor qLS(q) plotted versus ql p . Data are for semiflexible chains based on the SCLM and BFM including 4 choices of the 
stiffness each. In (a)(b), the straight lines indicate the rod-like behavior at large q (slope = —1) the SAW behavior for flexible 
chains (slope = — 1 /v, with v = 0.588), and the Gaussian behavior (slope = —2). In (c), the formulas proposed by Kholodenko 
{Eqs. (1451) - (I47|) }. the Debye function {Eq. (1441) 1 for Gaussian chains, and qLS(q) —> tt for a rigid-rod [57| are also shown for 
comparison. 


behavior causes slow transients in many quantities, be¬ 
fore the asymptotic behavior of flexible chains is reached 
(see e.g. Fig. ED- In the absence of the excluded vol¬ 
ume effect, a single crossover occurs, from rigid-rods to 
Gaussian coils as implied by the Kratky-Porod model, 
while a double crossover occurs from rigid-rods to Gaus¬ 
sian coils and then to swelling coils due to the excluded 
volume interaction as predicted by the Flory-like argu¬ 
ments. We have verified the Kratky-Porod crossover scal¬ 
ing behavior for semiflexible RWs, semiflexible NRRWs, 
and for semiflexible SAWs when the excluded volume ef¬ 
fect is not yet important, otherwise the Flory predic¬ 
tion takes over for semiflexible SAWs. The flexibility of 
chains in our model is controlled by the bending poten¬ 
tial U b = £f,(l — cos 9). Our results of bond-bond orienta¬ 
tional correlation functions (cos 9(s)) (Fig. fT^l) show that 
the persistence lengths of semiflexible RWs, NRRWs, and 
SAWs are the same for a given bending energy e b based 
on the same lattice model. But, with a caveat: there is 
a problem of fitting the exponential decay to the data of 
(cos 9(s)) for not only semiflexible SAWs but also semi¬ 


flexible RWs and NRRWs based on the BFM for e b > 10 
(rather stiff chains) due to the fluctuations of bonds and 
the lattice artifacts as it was mentioned in Ref. [H2] . The 
structure factor describing the scattering from semiflex¬ 
ible linear polymer chain based on the SCLM provides 
an almost perfect match to the result based on the BFM 
when we adjust the bending energy e b such that the same 
persistence length i v results for both models. 

From our simulations the different crossovers to the 
asymptotic behavior of single chains based on the SCLM 
and BFM are observed and investigated. Similar effects 
have to be expected for real chemical systems as well. 
Thus coarse graining will require different mapping ra¬ 
tios for different coarse-grained models. However, the 
equilibration time may rise dramatically for simulating 
large and complex realistic polymer systems. A proper 
mapping onto a coarse-grained model where the number 
of degrees of freedom is reduced should help to speed up 
the simulations. Based on the BFM, the bond angles and 
bond lengths of polymers can be treated as dynamic de¬ 
grees of freedom depending on temperature. Thus, the 
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static structure of a polymer model on the coarse grained 
level could be tuned, when one introduces bond length 
and bond angle potentials, to mimic the structure of a 
chemically realistic model of a polymer which contains 
covalent chemical bonds, whose orientation is controlled 
by both bond angle and torsional potentials. In this pa¬ 
per we did not discuss the details of this mapping pro¬ 
cedure yet, but we hope that our work will be a useful 
input for this problem. However, it will also be inter¬ 
esting and important to understand the distributions of 
bond lengths and torsional angles. 

We hope that the present work will contribute to a bet¬ 
ter understanding of using the lattice models for studying 
complex polymer systems and for the development of a 
multi-scale coarse-graining approach based on the lattice 


models. 
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